International Journal of Mathematics and 
Computer Applications Research (UMCAR) 
ISSN(P): 2249-6955; ISSN(E): 2249-8060 
Vol. 4, Issue 5, Oct 2014, 29-40 
© TJPRC Pvt. Ltd. 




UNSTEADY MHD FLOW THROUGH POROUS MEDIUM PAST A 



MOVING VERTICAL PERMEABLE SURFACE IN THE PRESENCE OF HEAT 



SOURCE AND A CHEMICAL REACTION 



P. R. SHARMA & MANISHA SHARMA 



Department of Mathematics, University of Rajasthan, Jaipur, Rajasthan, India 



ABSTRACT 



An analysis is made to present unsteady flow, heat and mass transfer of an electrically conducting, viscous and 
heat generating / absorbing fluid through porous medium past a uniformly moving vertical permeable surface in the 
presence of a magnetic field and heat source considering the first order homogeneous chemical reaction and energy loss 
due to viscous and Joulean heat dissipations. The effects of various physical parameters on velocity, temperature and 
concentration profiles are shown through graphs and coefficient of skin - friction, Nusselt number and Sherwood Number 
are presented through tables. 
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INTRODUCTION 

Flow through porous media is encountered in many branches of engineering and sciences e.g. ground water 
hydrology, reservoir engineering, solid sciences, soil mechanics and chemical engineering etc. Recently attempts have 
been made to study the effects of magnetic field on an electrically conducting fluid for its various applications in MHD 
power generators, intercontinental ballistic missiles, confinement of plasma in nuclear fusion etc. Soundalgekar considered 
effects of mass transfer and free convection currents on the flow past an impulsively started vertical plate[20]. Sharma and 
Mathur studied the steady laminar free convection flow of an electrically conducting fluid along a porous hot vertical plate 
in the presence of heat source / sink[13]. Anjali devi and Kandasamy observed the effects of chemical reaction, heat and 
mass transfer on a laminar flow along a semi-infinite horizontal plate[3]. Unsteady two dimensional laminar flow of a 
viscous incompressible electrically conducting fluid in the vicinity of a semi-infinite vertical porous moving plate was 
discussed by Kim [11]. Das et al. reported heat transfer effects on the flow of a viscous, incompressible and electrically 
conducting fluid between two stretched /squeezed horizontal porous plates subject to uniform injection at upper plate in the 
presence of a transverse magnetic field[7]. The effects of mass transfer on unsteady MHD flow and heat transfer past an 
infinite vertical porous moving plate was analyzed by Sharma and Mishra [14]. Israel - Cookey et al. observed the 
influence of viscous dissipation and radiation on unsteady MHD free convection flow past an infinite vertical plate in an 
optically thin environment with time dependent suction[10]. Hydromagnetic heat and mass transfer in the flow of a viscous 
incompressible fluid past an infinite vertical porous plate under oscillatory suction velocity normal to the plate was 
discussed by Singh et al. [19]. Afify studied the effects of a chemical reaction on the free convective flow and mass 
transfer of a viscous, incompressible and electrically conducting fluid over a stretching surface[l]. Chen discussed the 
combined heat and mass transfer effect on an electrically conducting fluid in MHD natural convection adjacent to a vertical 
surface taking into account the effects of Ohmic heating and viscous dissipation[5]. Sharma and Mishra presented unsteady 
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flow and heat transfer along a porous vertical surface bounded by porous medium[18]. Sharma and Gupta discussed 
unsteady flow and heat transfer along a hot vertical porous plate in the presence of periodic suction and heat source[17]. 
Das observed the effects of constant suction and injection on MHD three dimensional Couette flow and heat transfer 
through porous medium[6]. Mass transfer effect on unsteady mixed convective flow and heat transfer along an infinite 
vertical plate bounded with porous medium was investigated by Sharma and Katta [16]. Gangadhar presented radiation and 
viscous dissipation effects on chemically reacting MHD boundary layer flow of heat and mass transfer through a porous 
vertical flat plate[8]. Sharma et al. studied mass transfer with chemical reaction in MHD convective flow along a vertical 
stretching sheet[15]. A sinusoidal fluid injection / suction on MHD three dimensional Couette flow through a porous 
medium in the presence of thermal radiation was discussed by Ahmed and Kalita [2]. 

The objective of the present paper is to investigate the effects of heat and mass transfer of a unsteady MHD flow 
of an electrically conducting, viscous and heat generating / absorbing fluid through porous medium past a uniformly 
moving vertical permeable surface in the presence of a heat source considering the first order homogeneous chemical 
reaction and energy loss due to the viscous and Joulean heat dissipations. 

MATHEMATICAL ANALYSIS 

Heat and mass transfer on unsteady, laminar and two dimensional incompressible flow of an electrically 
conducting viscous fluid through porous medium, whose surface is maintained at a uniform temperature and uniform 
species concentration near the wall, are considered. The surface is assumed to be infinitely long i.e. the dependent variables 
do not dependent on vertical or axial coordinate. The fluid flows over a continuously moving vertical permeable surface 
with surface suction and heat generation / absorption. A uniform transverse magnetic field is applied normal to the 
direction of the flow. 

The magnetic Reynolds number is assumed to be small so that the induced magnetic field is neglected. The 
external electric field is zero and the electric field due to polarization of charges is negligible. Moreover a chemically 
reactive species is emitted from the moving vertical surface in hydrodynamic flow field. It diffuses into the fluid where it 
undergoes a simple isothermal, homogeneous chemical reaction. The reaction is assumed to take place entirely in the 
stream. Under these assumptions along with the Boussinesq approximation, the governing equations of unsteady flow, heat 
transfer and mass transfer in the presence of viscous and Joulean heat dissipation are given by 

dv 

— - = 0 =^> v is independent of y , 



du * du* d 2 u • • 

— — +gp(T -T m ) + gfi(Q -C x ) 
at dy dy 



aB 0 2 v 
— + 

P K 0 J 
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dt* dy' 
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where x* is taken along the motion of the surface in upwards direction and y is normal to it, u andv* are the 

velocity components along the directions of x and y , respectively, V the kinematic viscosity, g the acceleration due 

to gravity, (5 the coefficient of thermal expansion, f3 the coefficient of concentration expansion, T the temperature of 

the fluid, the temperature of the fluid far away from the wall, C the mass concentration, the mass concentration 

far away from the wall, (J the electrical conductivity of the fluid, B Q the magnetic field intensity, p is the fluid density, 

K Q the permeability of the porous medium, K the thermal conductivity, C p the specific heat at constant pressure, S the 

heat source / sink parameter, jU the coefficient of viscosity, D the mass diffusion coefficient and ?C the reaction rate 

term. Following Aris [4], reaction rate term X represents the reaction kinetics of the system whose overall reaction is 
described by the power-law model given by 

A m =-K\C*-CJ n , (5) 

where K is the reaction rate constant and n is the order of the reaction. Hence eq. (4) is reduced to 



*3C* 8 2 C* 
oy dy 



v— = D^^-K\C*-CJ, (6) 



The boundary conditions are 

y*=0 :u=U w , v=-V w (l + ee iM ), T* = T W , C*=C W ; 

(7) 

where U w (a constant), V w > 0 and T w are the surface velocity, mean suction velocity, temperature at the wall, 
respectively and C w is concentration at the wall. Solution of the equation (1) subject to equation (7) is 

v=-V w (l + £e i<or ). 

Introducing the following non-dimensional quantities 

i 

r -t,2 at 



uv- 



pcy; vi ' v 2 K 



Ec = — , Sc = — , ^ = ££. 

C P {T V -TJ D Vl 

(8) 

in the equations (2), (3), and (6) we get 
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d 2 u du du 
dy 2 dy dt 
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dd 
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K 



u — 0, 



(9) 
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Pr Pr — + SPrd + Ec Pr 
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+ M 2 EcVxu = 0, 



d 2 C BC 8C 

^-Sc — -Sc — -KScC = 0. 

dy dy dt 



(10) 



(ID 



where Gr the Grashof number, Gm the modified Grashof number, Pr the Prandtl nmber, Sc the Schmidt 
number, K chemical reaction parameter, M the Hartmann number, S heat source / sink parameter, K p permeability 

parameter, Ec suction Eckert number and other physical quantities have their usual meanings. 

The non-dimensional boundary conditions are given by : 



3; = 0 : u = l, v = -(\ + £e' m ),9 = \,C = \; 
3;^oo :M ^0, 0^0, C^O. 



(12) 



Equations (9) to (11) are coupled, non linear partial differential equations and their closed form solution cannot be 
determined. In the neighbourhood of the wall, velocity, temperature and concentration of the fluid can be expressed in the 
given form: 

u(y, t) - u 0 00 + se im u x ( y) + o(s 2 ) + . . 
0(yA) = e o (y) + ee«"e i (y) + o(e 2 )+... 
C(y, t) = C 0 (v) + ee ia 'C x (y) + o{s 2 )+.... 

substituting (13) in to equations (9) to (11) and equating the harmonic and non-harmonic terms and neglecting the 
higher order terms of 0(s) 2 , we obtain six equations for U Q9 U 19 # 0 , 0 19 C Q9 C x which are still coupled and non-linear 
whose exact solutions cannot be determined. Further, for incompressible fluid flows, the Eckert number is very small, 
therefore U Q , U x , 6 Q , 6 X , C 0 , C { can be expanded in the power of Ec as given below 



F{y)=F 0 {y)+EcF l {y)+O{Ec 1 ), 



where F stands for u Q9 U 19 0 Q9 0 19 C o orC v 



(14) 



Using eq. (14) into the expression of U Q9 U 19 0 Q9 0 l9 C Q9 C x and equating the terms free from Ec and coefficient 
of Ec , we get the following systems of ordinary differential equations: 
Zeroth order equations 



M 2 + 



-Gr0 oo -GmC oo , 



(15) 
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W 01 ^01 



f M 2 + — 



v 



a: 



»oi = -Gr6 0l - GmC 0l 



p J 



61 + Pr 0' Oi + S Pr 0 O1 = - Pr(u m f-M 2 Pr 2 u 2 m , 



C'^ + ScC' 00 -KScC 00 =0, 



C" { + Sc C' Ql -K Sc C 0] = 0. 



The corresponding boundary conditions are 

y = 0: u 00 =l, u 01 = 0, 0 OO = 1, 0 Q1 = 0, C 00 = 1, C 01 = 0; 

y = 0: W(X) =1 3 i/ 01 =0, 0^=1, 0 O1 = 0, C 00 =LC 01 =0;T 

v -> qo: Woo -> 0, z/ 01 -> 0, 6> 00 -> 0, 0 O1 -> 0, C 00 -> 0, C 01 -> 0. 

First order equations 



u l0 +u w 



10) + 



M 2 +- 



K. 



p J 



*10 "00 



"oo- Gr6> io- GwC io 5 



1 



10) + 



1 



e; o +?ra; o -?r(s-w))0 lo =-?r0' oo , 

6>n + Prd' n -Pr(S-io))0 n =-Pr^ 1 -2Pr»>[ 0 -2M 2 Pra^, 

c;; +&c; 0 -&(^+^)c 10 = -&c;, 

C" +SCC,', -&(wa+i:X^, = -&Coi. 

The corresponding boundary conditions are 

y = 0 : m 10 = 0, u n = 0, 0 W = 0, 0 n = 0,C 10 = 0, C n = 0; 

y -> «: a 10 -» 0, « n -> 0, <? 10 0, <? n -> 0, C 10 -> 0, C n -> 0. 
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Solving equations (15) to (20) under the boundary conditions (21) and the equations (22) to (27) under the 
boundary conditions (28), the expressions of W 00 , W 01 , 0 m , 0 Ql , C 00 ,C 01 , U l0 , W n , & l0 , 0 n ,C w and C n are known. 

Finally in view of eq. (13), the expressions of u(y, t) , 6 (y,t) and C(y,t) are known and given by 

u(y, t) = [(B 2 e Uy - B 3 e hy - B 4 e hy ) + Ec{B n e hy + B X3 e hy + B l4 e~ 2hy + B l5 e~ 2,iy + B l6 e~ 2hy + B l7 e~ Hiy 

+B n e-" 2y + B 19 e-" iy )] + £e im \{B 2x e hy + B 22 e~' 2y + B 23 e hy + B 24 e hy + B 2i e hy + B 26 e~ hy ) 

+Ec(B 44 e~ hy + B 45 e- hy + B 4(s e~ hy + B 41 e~ 2hy + B 48 e" 2, ' v + B 49 e 2hy + B 50 e" fl ' y + B,^ + B 52 e~" iy + B 53 e~"< y 

+B 54 e~" sy + B 55 e~ H6y + B i( e~ Hiy + B 51 e~ H * y + B 5S e~"' y + 5 59 e" Hl(,v + 5 60 e~ H " v + B 6l e~" liy ) ...(29) 

0( y ,t) = [e"' 3y + J Ec(B 5 e-' 3y + B^ + B^ + B^- 2hy + B 9 e~ Hiy + B,^ + B n e H ' y )] 

+£e ,m [(B 20 e~ hy - B 2Q e~ hy ) + Ec{B 21 e~ hy + B n e hy + B 29 e 2hy + B^e 2hy + B 31 e~ 2 ' 33, + 5 32 e~ Hl3, 

+B 33 e"^ + B M e H * + B 35 e~ H " y + B 36 c H > y + B 31 e H * y + B 3% e^ y + B 39 e~ H ' y + B 40 e^ y + B 4l c H « y 

+B 42 e~ H " y + B 43 e~ H " y )], ...(30) 

C(y,t) = e~ hy +ee im \-B x e hy + B x e~ hy ], ...(31) 

Skin - Friction Coefficient 

Knowing the velocity field, the skin - friction coefficient at the surface is given by 



T w _ ( du. _ du* 

PUJ W dy dy 



[B 62 +EcB 6i ] + se ta! [B 64 +EcB 6S ]. 



(32) 



Nusselt Number 

Knowing the temperature field, the rate of heat transfer In terms of Nusselt number at the surface is given by 
Nu = ?Z • q 

= [-l i + EcB 66 ] + £ e""[B 61+ EcB 6S ]. 
Sherwood Number 

Knowing the concentration field, the rate of mass transfer in terms of Sherwood number at the surface is given by 

= -I { +se B 69 . 

(34) 

Here, B { to B 69 , 1 { to I 6 and H x to H u are constants and their expressions are not given here for sake of brevity. 
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RESULTS AND DISCUSSIONS 

To study the physical situation of the problem, we have computed the numerical values of the velocity, 
temperature, concentration, skin -friction coefficient, rate of heat transfer in terms of Nusselt umber and rate of mass 
transfer in terms of Sherwood number for different values of the physical parameters. For velocity, temperature, 
concentration profiles skin - friction coefficient, Nusselt number and Sherwood number, parameter £ is valued as .25, 

CO is valued as 5, and OX is valued as 71 1 3 . 

It is observed from Figure 1 that fluid velocity decreases due to increase in the Hartmann number, while it 
increases due to increase in the Eckert number or heat source parameter. It is seen from Figure 2 that the fluid velocity 
decreases with increase of chemical reaction parameter, Schmidt number or Prandtl number. It is observed from Figure 3 
that fluid velocity increases due to increase in permeability parameter, Grashof number or modified Grashof number. 
Figure 4 shows that fluid temperature increases due to increase in heat source parameter, while it decreases due to increase 
in chemical reaction parameter or Schmidt number. Figure 5 illustrates that fluid temperature increases due to increase in 

Grashof number, modified Grashof number or permeability parameter where 1 < K < 2 . It is seen from Figure 6 that 

fluid temperature increases due to increase in Eckert number, while it decreases due to increase in Prandtl number or 
Hartmann number. It is seen from Figure 7 that mass concentration decreases due to increase in chemical reaction 
parameter or Schmidt number. 

It is observed from Table 1 that skin -friction coefficient increases due to increase in Grashof number, modified 
Grashof number, heat source parameter, Eckert number or permeability parameter, while it decreases due to increase in 
Hartmann number, Chemical reaction parameter, Schmidt number or Prandtl number. Nusselt number increases due to 
increase in Prandtl number, chemical reaction parameter or Schmidt number, while it decreases due to increase in Grashof 
number, modified Grashof number, Hartmann number, heat source parameter, permeability parameter or Eckert number. 
It is also seen from the Table 1 that Sherwood number increases due to increase in chemical reaction parameter or Schmidt 
number. 

CONCLUSIONS 

In view of the graphs and tables, the following conclusions are made: 

• Magnetic field produces a drag force called Lorentz force which causes reduction in the fluid velocity. 

• Buoyancy effect due to Grashof number and modified number sets a convection current which accelerates the 
velocity of the fluid particles. 

• Fluid velocity decreases with an increase in chemical reaction parameter, Schmidt number or Prandtl number. 

• The velocity increases due to increase in heat source parameter, Eckert number or permeability parameter. 

• Fluid temperature increases with an increase in heat source parameter, Grashof number, modified Grashof 
number, permeability parameter or Eckert number. 

• An increase in Hartmann number, chemical reaction parameter, Schmidt number or Prandtl number leads to 
decrease in fluid temperature. 
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• An increase in Schmidt number or chemical reaction parameter leads to decrease in mass concentration. 

• Grash of number, modified Grash of number, heat source parameter, Eckert number or permeability parameter 
enhance the coefficient of skin - friction whereas the reverse effect is observed for Hartmann number, chemical 
reaction parameter, Schmidt number or Prandtl number. 

• Nusselt number increases with an increase in Prandtl number, chemical reaction parameter or Schmidt number 
whereas it decreases with an increase in Grashof number, modified Grashof number, heat source parameter, 
permeability parameter, Eckert number or Hartmann number. 

• Chemical reaction parameter or Schmidt number enhance the Sherwood number. 



Table 1: Numerical Values of Skin-Friction Coefficient, Nusselt Number and 

Sherwood Number at the Surface for Various Values of Physical Parameters 
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Figure 2: Velocity profiles for different values of K, Sc and Pr, when M=l, 
S=0.1, ffl -5, K p =l, Ec=0.01, Gr=2, Gm=l, £ =0.25 and 0*=JT/3 
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Figure 3: Velocity Profiles for Different Values of K p , Gr and Gm, when S=0.1, 
® =5, Sc=0.22, Ec=0.01, M=l, K=2, Pr=0.71, £ =0.25 and &t =xl 3 
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Figure 4: Temperature Profiles for Different Values of S, K and Sc, when ® z =5, 
Ec=0.01, M=l, Gr=5, Gm=4, K p =l, Pr=0.71, £" =0.25 and ^=ff/3. 
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Figure 5: Temperature Profiles for different values of Gr, Gm and Kp, when S= 0.1, 
& =5, K=2, Ec=0.1, m=l, Sc= 0.22, Pr=0.71, £ =0.25, ^ = ff/ 3 - 




Figure 6: Temperature Profiles for Different Values of M, Ec and Pr, when S=0.1, 
® =5, K=2, Gr=5, Gm=4, Sc=0. 22, K p =l, ^=0.25 a nd &*=xf3 




Figure 7: Concentration Profiles for Different Values of K and Sc, when S=0. 1, 

0 =5, Ec=0.01, Gr=2, Gm=l, M=l, K p =l, Pr=0.71, £ =0.25 and M = ff/ 3 - 
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